Hole— defect chaos in the one— dimensional complex Ginzburg— Landau equation 
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We study the spatiotemporally chaotic dynamics of holes and defects in the ID complex Ginzburg- 
Landau equation (CGLE). We focus particularly on the self-disordering dynamics of holes and on 
the variation in defect profiles. By enforcing identical defect profiles and/or smooth plane wave 
backgrounds, we are able to sensitively probe the causes of the spatiotemporal chaos. We show that 
the coupling of the holes to a self-disordered background is the dominant mechanism. We analyze a 
lattice model for the ID CGLE, incorporating this self-disordering. Despite its simplicity, we show 
that the model retains the essential spatiotemporally chaotic behavior of the full CGLE. 
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I. INTRODUCTION 

The formation of local structures and the occurrence 
of spatiotemporal chaos are among the most striking 
features of pattern forming systems. The complex 
Ginzburg-Landau equation (CGLE) 

d t A = A + (1 + i Cl )d 2 x A - (1 - ic 3 )\A\ 2 A (1) 

provides a particularly rich example of these phenomena. 
The CGLE is the amplitude equation describing pattern 
formation near a Hopf bifurcation 0, , and exhibits an 
extremely wide range of behaviors as a function of c\ and 

ca Quassia. 

Dejects and holes are local structures that play a cru- 
cial role in the intermediate regime between laminar 
states (small ci, C3) and hard chaos (large ci, C3). Iso- 
lated defects occur when A goes through zero, where the 
complex phase ip := arg(A) is no longer defined. Ho- 
moclinic holes are localized propagating "phase twists" , 
which are linearly unstable. As illustrated in Fig.^ holes 
and defects are intimately connected: defects can give 
rise to "holes", which may then evolve to generate de- 
fects, from which further holes can be born, sometimes 
generating self-sustaining patterns. For more details see 
Refs. H0. 

The aim of our paper is to understand and model the 
spatiotemporally chaotic hole-defect behavior of the ID 
CGLE, built on the local interactions and dynamics of 
the holes and defects. Given the strength of the initial 
phase twist that generates a hole, and the wavenumbcr 
of the state into which it propagates, the hole lifetime r 
turns out to be the key feature. Surprisingly, the initial 
phase twist and invaded state play very different roles. 
For hole-defect chaos, we will show that the defect pro- 
files, which constitute the phase twist initial condition 
for the resulting daughter holes, show rather little scat- 
ter for fixed ci and C3. Changes in c\ and C3, however, are 
encoded in changes in the defect profiles, and thus lead 
to changes in the typical lifetimes of the daughter holes. 
We then demonstrate that the chaos does not result from 
variations in defect profiles. It rather follows from the 
sensitivity of the holes to the states they invade, since the 
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FIG. 1: Illustration of the main phenomenology of hole-defect 
chaos (after Refs. 0,0,0). (a-b) Space-time grey-scale plots 
showing the invasion of a plane wave state by hole-defect 
chaos: (a) \A\ (dark: AsiO) and (b) q := d x ip (light: q > 0, 
grey: g~0, dark: q<0). The propagating objects are incoher- 
ent holes, which dynamically connect the defects (the black 
dots in (a)). Parameter values are c\ = 0.6, C3 = 1.4 , with 
an initial condition given by Eq. with 7 = 1, q e x = —0.03. 
The non-zero q ex breaks the left -right symmetry and results 
in the differing periods of the left and right moving edge holes, 
(c) Close-up of I .A I and (d) close— up of the complex phase tjj, 
showing in detail how a hole generates a phase defect that in 
turn generates two daughter holes. 



passage of each hole disorders the background wavenum- 
ber profile leading to disordered background states. It is 
the self-disordering action of the holes that is primarily 
responsible for the spatiotemporal chaos. 

With these insights, we can then construct a simpli- 
fied lattice model for hole-defect chaos, which both re- 
produces the correct qualitative behavior as C\ and C3 are 
varied and which captures the correct mechanism (prop- 
agating, self-disordering holes). Our initial findings on 
this subject can be found in Ref. 0, where we introduced 
the concept of self-disordering, and outlined a simplified 
lattice model. However in this paper, we investigate the 
subject in considerably greater depth, and, in particular, 
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provide much more conclusive evidence for the correct- 
nesss of the self-disordering hypothesis. 

We now briefly summarize the structure of the pa- 
per. Topics discussed already in earlier work 0, Q are 
dealt with rather briefly. We start in Section II by de- 
scribing hole-defect dynamics on a local scale. In Sec- 
tion III, we then use this knowledge to investigate disor- 
dered hole-defect dynamics, and conclusively show that 
it is the coupling of the holes to a self-disordered back- 
ground which is the dominant mechanism for spatiotcm- 
poral chaos. This concept is then illustrated by a mini- 
mal lattice model for hole-defect dynamics in Section IV, 
before we draw our conclusions in Section V. 



II. HOLE-DEFECT DYNAMICS 

We begin by studying the hole lifetime r as a function 
of the initial conditions (Fig. [2J. This study motivates 
the central question of this paper: how does r depend on 
the initial conditions and on the external wavenumber, 
and which of these dependencies is most important for 
spatiotemporal chaos? We then study general properties 
of defect profiles, and demonstrate that in hole-defect 
chaos the profiles of defects show rather little scatter. 
We also show how the lifetimes of "daughter" holes born 
from a typical defect vary with c\ and C3. Taken to- 
gether, the data presented here forms direct evidence for 
the heuristic picture of hole-defect dynamics developed 
in Refs. %%. 

A. Incoherent homoclons 

In full dynamic states of the CGLE, one does not ob- 
serve the unstable coherent homoclinic holes unless one 
fine-tunes the initial conditions (see below). Instead 
evolving incoherent holes which either decay or start to 
grow towards defects occur 0, 0- Let us consider the 
short-time evolution of an isolated incoherent hole prop- 
agating into a regular plane wave state. Holes can be 
seeded from initial conditions like Q: 

A = cxp{i[q EX x + (n/2) tanh^x)]) . (2) 

The two essential parameters 7 and q ex represent, re- 
spectively, the initial conditions from which the incoher- 
ent hole is born and the background wavenumber of the 
state into which the hole propagates. In this context, a 
single parameter 7 is sufficient to scan through different 
initial conditions, since the coherent holes have just one 
unstable mode J5J. 

A detailed contour plot of the lifetime of an initial inco- 
herent hole as a function of 7 and q ex is shown in Fig. 
These results were obtained using a semi-implicit nu- 
merical integration of the CGLE, with space and time 
increments dx — 0.25 and dt = 0.01. As expected, three 
possibilities can arise for the time evolution of the ini- 
tial peak: evolution towards a defect (upper right part of 
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FIG. 2: Contour plot showing the lifetime r of an initial in- 
coherent homoclon before a defect is generated. The initial 
condition is given by Eq.(|5J, and the lifetime is plotted as a 
function of 7 and q ex . Note that the lifetime diverges as q ex or 
7 are reduced. In the left-bottom corner of the diagram, the 
incoherent homoclon decays and no defects are formed. The 
inset shows a sketch of the phase space around the homoclon 
saddle (after Ref. 7]), where the manifold F represents the 
family of peaked initial conditions of the form 

Fig. I2J), decay (lower left part of Fig. [2J, or evolution ar- 
bitrary close to a coherent homoclinic hole (the boundary 
between these two regions). These possibilities, together 
with an illustrative sketch of the phase space, are shown 
in the inset of Fig. [21 The rather simple and monotonic 
behavior of r with q ex and 7 is somewhat of a surprise, 
and this reinforces our simple phase space picture; no 
other solutions seem to be relevant in this region of phase 
space. 

Since homoclons are neither sinks nor sources, Fig. [21 
can be interpreted as follows: for a right -moving homo- 
clon an incoming wave with positive wavenumber tends 
to push the homoclon more quickly towards a defect; 
previously we have referred to this as "winding up" of 
the homoclon. Similarly, an incoming wave with nega- 
tive wavenumber "winds down" a right-moving homo- 
clon, possibly even preventing the formation of a defect 

a 

B. Defects 

We now study the defect profiles themselves in more 
detail. In Fig.[3K we show complex plane plots of Re(a) 
vs. Im(A) just before, close to, and just after a defect. As 
can be seen, there is no singular behavior whatsoever: the 
real and imaginary parts are smooth functions of x and 
t, even at the time of defect formation. However, when 
transforming to polar coordinates, a singularity manifests 
itself at the defect, where \A\ — > 0. This can also be 
seen from the g-profiles shown in Fig. |2j3-d. In fact, it is 
straightforward to show that the maximum value of the 



3 




¥ v ¥ 

~ (e) ~ (f) 

-1 Re(A) 1 -1 Re(A) 1 -1 Re(A) 1 



1 <S> 



FIG. 3: (a) Plots of the real and imaginary part of A around 
a defect, just before (thin line), close to (medium line) and 
after (thick line) the defect has occurred; the time difference 
between successive profiles is 0.01. Note that in this represen- 
tation the dynamics looks completely smooth, (b-d) g-profilcs 
for these three cases (identical vertical scale in each case); just 
before the defect occurs, a large positive spike develops in the 
profile, while after the defect this becomes a large negative 
spike, (e-g) Statistics of defects obtained by overlaying 10 3 
defect profiles of spatial extension (width) 20, centered around 
the defect position. An arbitrary phase factor has been di- 
vided out by requiring that Re(d x (A) | defect) =0. All data was 
collected in a system of size 500, after a transient of 500. The 
coefficients a and c 3 are: (e) 0.6, 1.4, (f) 1.4, 1.4, (g) 3.0,3.0. 



local phase gradient q m diverges as (At) 1 at a defect , 
where At is the time before defect nucleation (see Fig. 3b 
of Ref. 0) 0. 

In Fig. EI e-g), we overlay complex plane plots of A 
around 10 3 defects obtained from numerical simulations 
of the CGLE in the chaotic regime. Surprisingly (see 
Fig. , defect profiles of the interior chaotic states are 
dominated by a single profile in the hole-defect regime, 
similar to fixing 7 in Fig. [21 This provides a strong indi- 
cation that hole-defect chaos does not come from scat- 
ter in the defect profiles. For large enough ci and C3, 
where holes no longer play a role, and where hard defect 
chaos sets in the profiles show a much larger scatter 
(Fig. Kg)). 
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FIG. 4: (a) Example of the divergence of a near-coherent 
right-moving hole, showing also the definition of times ti and 
<2. (b) Times t\ (diamonds) and ti (triangles) as a function 
of C3 for ci = 0.6. The inset shows that the curves for ti and 
ti actually cross for sufficiently large C3. 



Hence the defect profile acts as an initial condition for 
its daughter holes, as can also clearly be seen in Fig. 

We now examine the fate of these daughter holes in 
the well-defined case where the initial defect is generated 
from the divergence of a right-moving, near-coherent ho- 
moclon in a q ex — background state (see Fig. 0Jt). We 
then define t\ and ti as the lifetimes of the resulting 
daughter holes. When a daughter hole does not grow out 
to form a defect, its lifetime diverges. In Fig.^J) we plot 
t\ and ti for c\ = 0.6 as a function of C3. The initial hole 
that formed the first defect has a lifetime of at least 60, 
and we have checked that a further increase of this time 
does not change t\ and ti appreciably. When both t\ and 
ti are infinite, no defect sustaining states can be formed, 
and the final state of the CGLE is in general a simple 
plane wave. When only t\ is finite, isolated zigzag states 
are formed; such states have been discussed in Ref. 
and we will see some examples below. When both t± and 
ti are finite, and of comparable value, more disordered 
states occur. We will later use this data on t\ and ti to 
calibrate our minimal lattice model for spatiotemporal 
chaos. 

Hence, we see that changes in Ci and C3 not only lead 
to changes in the defect profiles, but also modify the life- 
times of the resulting daughter holes. However, for fixed 
ci and C3, we have seen that the defect profiles, which act 
as initial conditions for the daughter holes, show rather 
little scatter, at least in the hole-defect regime. In the 
next section, we will build on this knowledge to unravel 
the causes of hole-defect spatiotemporal chaos. 



C. Defect — > holes 

Suppose a hole has evolved to a defect; what dynamics 
occurs after this defect has formed? As Fig. [3{ b-d) shows, 
defects generate a negative and positive phase-gradient 
peak in close proximity. The negative (positive) phase 
gradient peak generates a left (right) moving hole, and 
analogous to what we described in Fig. |3 the lifetimes 
of these holes depend on the initial peak and on q ex . 



III. MECHANISM OF HOLE-DEFECT CHAOS 

In this paper and in earlier work 0, we have argued 
that the principal cause of the spatiotemporally chaotic 
behavior in the ID CGLE is the movement of holes 
through a self-disordered background. Clearly, as we can 
see from Fig.0 a disordered background wavenumber q ex 
will give rise to varying hole lifetimes and thus to disor- 
dered hole-defect dynamics. In this section, we explicitly 
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FIG. 5: (a) Space-time plot of q in the ordinary CGLE for 
ci = 0.6, C3 = 1.4. (b) Log-greyscale plot of growth of per- 
turbations. At t = 166.66, x = 500, one gridpoint was altered 
by 10~ 6 . (c) Space-time plot of q in model (i), the CGLE 
with fixed defect profiles, (d) Log-greyscale plot of growth of 
perturbations for the fixed defect model (i). 



demonstrate the correctness of this mechanism by modi- 
fying the CGLE dynamics in two ways: 
Model (i): Fixed defect profile. Whenever a defect occurs, 
this defect is replaced by a standardized defect profile 
(obtained from an edge defect). Here the dynamics will 
be chaotic, showing the irrelevance of the scatter in defect 
profiles. 

Model (ii): Background between holes — > plane wave with 
q = 0. At each timestep, the background between any 
two holes is replaced by a plane wave with wavenumber 
zero. Here no chaos will occur, illustrating the crucial 
importance of the self-disordered background. 

In case (i), the size of the replaced defect profile was 
five centered around the defect; in case (ii) the back- 
ground was defined to be all regions where \A\ > 0.95. 
Our results are substantially independent of the exact de- 
fect size or cutoff value. In both cases, it is crucial to en- 
sure that no jumps in the phase occur at the edges of the 
replaced regions. This can be achieved by phase match- 
ing the replaced region (either defect profile or plane 
wave) at the left boundary, while the state to the right 
of the replaced region is multiplied by a phase factor to 
enforce phase-continuity at the right edge. We take open 
boundary conditions (i.e. d x A = 0) and only study the 
behavior far away from these boundaries. 

In Fig. [S] we show an example of the dynamics and 
spreading of a localized perturbation for the full CGLE 
(Fig.GJi-b) and for the "fixed defect" model (i) (Fig.Efc- 
d), both for c\ = 0.6, c 3 = 1.4. For both models, we took 
as an initial condition a defect rich state, which after 
a few timesteps shows the typical hole-defect dynamics. 




FIG. 6: Dynamical states in model (ii), i.e., the CGLE where 
the background between holes is replaced by a q = plane 
wave, (a) For C3 = 1.4 and c\ — 0.6 only a few isolated zigzags 
occur, (b) When ci is increased to 0.7, more zigzags occur, 
but there is no chaos, (c-d) For ci = 0.8, a disordered tran- 
sient occurs (c) that eventually freezes into a quasiperiodic 
zigzag state without disorder (d). 



At t ss 167 we applied a local perturbation of strength 
10~ 6 to the middle gridpoint (corresponding to x = 500) 
and followed the evolution of both the perturbed and 
the unperturbed systems in order to follow the spread- 
ing of perturbations. For the full CGLE (Fig. [5^-b) , the 
perturbation spreads along with the propagation of the 
holes. We note that the initial growth of the perturba- 
tion manifests itself in slight "shifts" of the spatial and 
temporal positions of the defects. In particular when two 
holes collide, a strong amplification of the perturbations 
is observed. 

We can now compare this with the above fixed defect 
profile model (i) (Fig. HJ, d). Clearly, the replacement of 
the defects does not destroy the chaotic behavior of the 
system, as confirmed by the spreading of a localized per- 
turbation (Fig. 01), which propagates in a similar fash- 
ion to the full CGLE (Fig. ^p). This strongly indicates 
that variation in the defect profiles is not contributing 
in a major way to the spatiotemporally chaotic behavior 
of the full CGLE. We should also point out one subtlety 
here: due to the discretization of time, the times at which 
defect profiles are replaced are also discrctized, and one 
may worry whether this destroys the chaotic properties of 
the model. However, we have performed simulations for 
a smaller timestep (dt = 0.001) and found no qualitative 
difference. As we will see, this issue of discretization will 
play a more important role in the lattice model discussed 
in section Hvl 

Turning now to model (ii), where laminar regions of the 
CGLE are replaced by q — plane waves, we see that the 
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FIG. 7: Grid model geometry showing the sites (dots) and 
hole propagation directions (arrows). The update rule is de- 
fined within a 2 x 2 cell, mapping (<j)i, <j> r ) — » ((/>{, <j>' r ). "Active 
sites" where > <j> t represent holes, while "inactive sites" 
where \4>\ < <t>t represent the diffusive background. When 
both sites are inactive the relevant dynamics is phase diffu- 
sion: (j)' r = D(j>i+(1—D)(j>r ■ The value of D is fixed at 0.05 and 
is not very important. When both sites are active, two holes 
collide and merge: <p' r — <j}[ — (<j>i+<j> r )/2. When one site is 
active but smaller than cj>d, we implement the evolution [Tof: 
4>'r = + A(c/>; — <j> n — g<f> T ) (we assume here that we have 
a right moving hole, the case of a left moving hole follows 
by symmetry). Here A sets the time scale and can be taken 
small (fixed at 0.1). When (j> > cj>d a defect occurs and two 
new holes, i.e., active sites, are generated; for details see text. 



disorder is destroyed. This is illustrated in Fig. |SJ where 
we show examples of model (ii) dynamics for C3 = 1.4 
and ci =0.6, 0.7, 0.8. Clearly chaos is suppressed for c\ — 
0.6 and 0.7 (Fig. EJi-b) , with zigzag type patterns being 
especially dominant. For C\ = 0.8 (Fig. [U), the initial 
dynamics do appear disordered, but after a transient the 
system evolves to a regular zigzag state (Fig. Eli). 

From the behavior of models (i) and (ii) we conclude 
that the self-disordered background is an essential ingre- 
dient for hole-defect chaos, while scatter in the defect 
profiles is not. 



IV. LATTICE MODEL 

To further justify and test our picture of self- 
disordered dynamics, we will now combine the various 
hole-defect properties with the left-right symmetry and 
local phase conservation to form a minimal model of 
hole-defect dynamics. The model reproduces regular 
edge states, spatiotemporal chaos and can be calibrated 
to give the correct behavior as a function of c\ and cj. 
An earlier version of the model was presented in Ref. Q • 
However, as will become clear, we have now modified and 
improved the model, and also for the first time made di- 
rect comparisons with the full CGLE. 

From our earlier analysis (see also Ref. 0), we see that 
the following hole-defect properties must be incorporated 
in the model: (i) Incoherent holes propagate cither left or 
right with essentially constant velocity, (ii) Their lifetime 
depends on c±, C3, and the wavenumber of the state into 
which they propagate. When the local phase gradient 
extremum diverges, a defect occurs. (Hi) Each defect, in 



turn, acts as an initial condition for a pair of incoherent 
holes. 

In our lattice model we discretize both space and time, 
and take a "staggered" type of update rule which is com- 
pletely specified by the dynamics of a 2 x 2 cell (see 
Fig. UJ. We put a single variable 4> on each site, where 
<j) corresponds to the phase difference (the integral over 
the phase-gradient q) across a cell, divided by 2ir. Local 
phase conservation is implemented by 4>i + <P' r — & + 4>r , 
where the primed (unprimed) variables refer to values af- 
ter (before) an update. Holes are represented by active 
sites where \<j>\ > </> t ; here <j> plays the role of the internal 
degree of freedom. Inactive sites are those with \<f>\ < (j> t , 
and they represent the background. The value of the cut- 
off (j)t is not very important as long as it is much smaller 
than the value of <f> for coherent holes, and is here fixed 
at 0.15. Without loss of generality we force holes with 
positive (negative) <j) to propagate only from <pi (<f> r ) to 

The details of the translation of these rules into the 
model can be found in the caption of Fig.[7|and in Ref. 0, 
with one exception. A "defect" is formed when fa > 
Here we have adopted two alternative schemes. In the 
simplest case (A) (studied before in Ref. 0), we take 
4>' r — 4>ad, and = 0d — 1 — 4> a d ■ Here we completely fix the 
new holes. The factor —1 reflects the change in the total 
winding number associated with the defects. Notice that 
the overall winding number does not change by exactly 
— 1. This is because <fii + <fi r is usually slightly different 
from <f>d- As we will discuss below, to avoid breaking 
this "phase conservation" we have also studied case (B), 
where we take <p' r = <f> a d, and = <pi + 4> r — 1 — <j) a d- Here 
some (small) scatter in the defect profiles is allowed, but 
the change in the overall winding number is now strictly 
-1. 

The model does contain a large number of parameters, 
g, <p n , 4>d and 4> a d- We will first discuss the role of g and 
the difference between the two defect rules (A) and (B). 

In order for the model to reproduce the correct lifetime 
dependence of edge-holes in hole-defect states, the cou- 
pling of the holes to their background, g, should be taken 
negative (although its precise value is unimportant). For 
g = the hole lifetime r becomes a constant, indepen- 
dent of the cj) of the state into which the holes propagate, 
and moreover, the dynamical states are regular Sierpin- 
sky gaskets (Fig. ^t). For g<0 both the appropriate r- 
divergence and disorder occur (Fig. Hp), illustrating the 
crucial importance of the coupling between the holes and 
the self-disordered background. 

It turns out that the model with defect rule (A) is not 
strictly chaotic; sufficiently small perturbations do not 
always chaotically spread. However, with defect rule (B) 
implemented, infinitesimal perturbations do spread (see 
Fig-Etid). To understand this difference consider the fate 
of a tiny, localized perturbation. Holes will sweep past 
and be influenced by this perturbation, but since time is 
discrete, a sufficiently small perturbation will not lead to 
a change in the time at which a hole evolves to a defect. 
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FIG. 8: Illustration of the necessary ingredients for disorder 
and chaos in the grid model. In all cases, A = 0.1, (f> a d = 0.75, 
4> n = 0.59, (f) d = 1.01 (this leads to ti = 10, t 2 = 12, which 
is the situation in the full CGLE for ci = 0.6 and C3 = 1.5). 
(a) Without coupling to the background, g = 0, the model 
with defect rule (A) leads to regular Sierpinsky gaskets, (b) 
When g — —3, the model generates disordered states, that 
are not strictly chaotic (see text), (c) When defect rule (B) 
is applied, also for g — —3, the dynamics is truly chaotic, as 
illustrated by the spreading of perturbations (d). 



We have found that after a number of holes have passed 
over such a perturbation, it can actually be absorbed, 
so that no chaotic amplification occurs. It is therefore 
the combination of the discreteness of time and the fixed 
defect profiles which do not, strictly speaking, lead to 
chaos. By lowering A, this problem is diminished, but 
this makes the model far less effective computationally. 
Alternatively, we have found that defect rule (B) also 
circumvents this problem; perturbations can now never 
be absorbed, due to the nature of the defect rule (B). In 
this case the defect profile is not entirely fixed, but its 
scatter is still rather small: for A = 0.1 a typical scatter is 
of the order of 5%, and this diminishes as A is decreased. 
Therefore we can conclude that, in the continuous time 
limit of the lattice model, the scatter of the defect profiles 
is not necessary to obtain chaos. In the remaining part 
of the paper, we will use model (B) only. 

The self-disordering can be very clearly observed in 
the minimal model, since its update rules unambiguously 
specify which sites are "background" and which are 'ac- 
tive" . Two snapshots of the evolution shown in Fig. 
are plotted in Fig. [5] These snapshots clearly demon- 
strate how, after sufficient time has passed, the "inactive" 
background between the holes has become completely 
disordered. 

The essential parameters determining the qualitative 
nature of the overall state are (j) n , <j>d and 4>ad- These 
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FIG. 9: Illustration of the self-disordering in two snapshots 
of the field <f> in our lattice model. Note that the scale is 
such that only the background is clearly visible; the peaks 
correspond to active, hole-like states. The data are taken 
from the runs shown in Fig.^, at (a) t = 60 and (b) t = 300. 
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TABLE I: Times t\ and t 2 as obtained in the full CGLE, and 
appropriate coefficients <f> n and (f>d that reproduce these times 
in our grid model. 



parameters determine the amount of phase winding in 
the core of the coherent holes with q ex = {<j) n ), and in 
the new holes generated by the defects (4>d, 4> a d)- When 
varying the CGLE coefficients, these parameters change 
too, leading to qualitatively different states. In particu- 
lar, they determine the times t\ and ti that we already 
studied for the full CGLE in section lITCl We found that 
when 4> n and 4>ad are both decreased, t\ and £2 roughly 
remain the same. We have therefore kept 4> a d = 0.65, and 
varied 4> n and 4>d to tune the values of t\ and £2- Notice 
that the symmetry (or asymmetry) of the defect profile 
depends on <pd — 1; a value of <f>d < 1 typically promotes 
zigzag patterns. We have determined the appropriate 
values of <p n and <pd for three concrete cases, tabulated in 
Table I. Notice that the parameters for C3 = 1.5 precisely 
correspond to those used in Fig. |S] 

As can be seen in Fig. I1UI the agreement between the 
simple model and the CGLE is satisfactory, although 
clearly the CGLE displays richer behavior. Note that in 
the full CGLE, small perturbations of the background 
wavenumber evolve in a non-trivial manner. For ex- 
ample, a nonzero average background wavenumber in- 
troduces a drift of the phase perturbations in the back- 
ground between the holes • Since this phase dynamics 
is much slower than the hole-defect dynamics, we have 
chosen to ignore it in the grid model, and this accounts 
for the difference between Fig. HOfa , and Fig. 110b . 

Finally, we emphasize that the grid model allows us to 
disentangle the mechanism of hole-defect chaos, by en- 
abling us to completely control the behavior of defects 
and the coupling between holes and the laminar back- 
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FIG. 10: Comparison between the dynamics of the full CGLE 
(a,c) and our grid model (b,d), where t\ and £2 are matched 
according to Table I, with c\ = 0.6, and (a-b) C3 = 1.25 or 
(c-d) c 3 = 1.4. 

ground. The grid model also has the advantage of being 
possibly the simplest model that captures the essence of 
the self-disordered hole-defect spatiotemporal chaos. We 
also emphasize that we have, for the first time, carried 
out a detailed comparison between the full CGLE and 
the grid model, both in our analysis of the spreading of 
perturbations and in the calibration of the grid model 
as a function of c\ and C3. Given the simplicity of the 
model, the agreement with the full CGLE is striking. 

V. CONCLUSION 

In conclusion, we have studied in depth the dynamics 
of local structures in the ID CGLE. We have presented 
strong evidence that the origin of the chaotic behavior 
in the ID CGLE lies in the self-disordering action of the 
holes, rather than in the scatter of the defect profiles. 
Using this insight, we have then developed a minimal 
lattice model for spatiotemporal chaos, which, despite 
its simplicity, reproduces the essential spatiotcmporally 
chaotic phenomenology of the full CGLE. 

How general are these results? We conjecture that 
there are two crucial properties needed for hole-defect 
type chaos: propagating saddle-like coherent structures 



(the holes) and a "conserved" field (the phase field). Of 
course, the phase is not strictly conserved here due to the 
occurrence of defects, but (f dx ip) mod 2tt is a conserved 
quantity. It is this conservation that is weakly broken in 
our grid model for defect rule (A), but is preserved for 
defect rule (B). Only the latter is truly chaotic. The con- 
servation is also the underlying reason why an evolving 
hole leaves an inhomogeneous and self-disordered trail 
behind. Without such a conserved field, there is no rea- 
son for "self-disordering" to occur, and the holes then 
typically will exhibit a fixed lifetime, leading to Sierpin- 
sky gasket type patterns as is often the case in reaction- 
diffusion models E2 • A related scenario appears to occur 
in the periodically forced CGLE EH . Conserved fields of 
the type described here may be expected more gener- 
ally for systems undergoing a Hopf bifurcation, and can 
therefore be expected to also occur in Ginzburg-Landau 
type equations including higher order terms, and also in 
experiments. We have argued in Ref. |14|, that saddle- 
type structures like the homoclons here, may be much 
more general. This leads us to believe that the type of 
dynamics described here is not an artefact of the pure 
CGLE, but could be far more widespread. 

Our work opens up the possibility for quantitative 
studies of hole-defect and homoclon dynamics, states 
which have recently been observed in various convection 
experiments [3 Eg. We hope that our simple picture 
will advance these experimental studies of space-time 
chaos into the quantitative realm. Local dynamics of the 
type studied here, such as the dependence of lifetime on 
initial conditions |14| , or measurements of quantities like 
the daughter hole lifetimes t\ and t% should be accessi- 
ble in experiment, thereby circumventing the difficulties 
normally associated with characterizing fully-developed 
chaotic states. 

Finally, we mention another commonly observed type 
of spatiotemporal chaos occurring in systems when a pe- 
riodic state undergoes a certain symmetry breaking bi- 
furcation E3 • Mathematically, such systems may be de- 
scribed by a complex Ginzburg-Landau equation, cou- 
pled to a phase field. Such models are sometimes referred 
to as A — <fi models 18]. Theoretically, the role of holes 
and defects has not yet been studied in great detail for 
these systems, but the main ingredients for hole-defect 
chaos of the type described here appear to be present. 
We hope that our work will encourage further studies in 
this area. 

M.H. acknowledges support from Stichting FOM and 
from The Royal Society. 
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